home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / ZBRENT.DEM < prev    next >
Text File  |  1991-04-29  |  705b  |  39 lines
  1. PROGRAM d9r7(input,output);
  2. (* driver for routine ZBRENT *)
  3. CONST
  4.    n=100;
  5.    nbmax=20;
  6.    x1=1.0;
  7.    x2=50.0;
  8. TYPE
  9.    glnbmax = ARRAY [1..nbmax] OF real;
  10. VAR
  11.    i,nb : integer;
  12.    root,tol : real;
  13.    xb1,xb2 : glnbmax;
  14.  
  15. (*$I MODFILE.PAS *)
  16. (*$I BESSJ0.PAS *)
  17.  
  18. FUNCTION fx(x: real): real;
  19. BEGIN
  20.    fx := bessj0(x)
  21. END;
  22.  
  23. (*$I ZBRAK.PAS *)
  24.  
  25. (*$I ZBRENT.PAS *)
  26.  
  27. BEGIN
  28.    nb := nbmax;
  29.    zbrak(x1,x2,n,xb1,xb2,nb);
  30.    writeln;
  31.    writeln('Roots of bessj0:');
  32.    writeln('x':18,'f(x)':15);
  33.    FOR i := 1 to nb DO BEGIN
  34.       tol := (1.0e-6)*(xb1[i]+xb2[i])/2.0;
  35.       root := zbrent(xb1[i],xb2[i],tol);
  36.       writeln('root ',i:2,'  ',root:12:6,fx(root):14:6)
  37.    END
  38. END.
  39.